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Abstract 

We discuss approximation ideas that can be used in parameter estimation 
and feedback control for Euler— Bernoulli models of elastic systems. Focusing 
on parameter estimation problems, we outline how one can obtain convergence 
results for cubic spline— based schemes for hybrid models involving an elastic 
cantilevered beam with tip mass and base acceleration. Sample numerical 
findings are also presented. 
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I. Introduction 


In this lecture we discuss some approximation techniques that may 
be used in algorithms for parameter estimation and/or feedback control 
in distributed models such as those arising in models typical of large 
flexible space structures. The focus of our recent efforts has been the 
development and analysis of computational algorithms, e.g., convergence 
analysis, numerical implementation (software development) and testing. 
While the ideas involved are also applicable to the computation of 
feedback controls, we restrict our discussions here to some of our 
efforts on techniques in the context of parameter estimation or 
"inverse” problems: given observations of a system, determine 
parameters in models which best describe structural/material properties 
manifested by the system in response to perturbations (loading, etc.). 

The importance of such problems is twofold: (i) parameter 
estimation can be viewed as a primary tool in on-orbit model development 
and analysis where one seeks to understand elastic/ viscoelastic material 
properties such as damping, stiffness, etc. and to detect changes in 
these due to aging, prolonged stress, etc.; (ii) parameter estimation is 
a precursor to and integral part of development of sophisticated 
feedback control laws (via feedback operators satisfying infinite 
dimensional Riccati equations involving functional parameters of the 
system) • 

Many of the structures of interest to aerospace engineers entail 
systems composed of composite materials in rather complex 
geometric/structural configurations. The need for methods to 
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investigate such variable structure distributed models has, in our 
opinion, been clearly established in a number of recent efforts 
including [1], [2], [3], [4]. A number of investigations of parameter 
estimation in models for elastic beams have involved approximation 
results (the Trotter-Kato theorem) from linear semigroup theory. In 
particular, problems for simple beams have been treated in this manner 
in [5], [6], [7], [8]. In [9] the Trotter-Kato ideas are employed to 
establish results for hybrid models similar to those introduced later in 
this presentation and which are important in the study of shuttle- 
deployed payloads. However, in some instances it is advantageous to use 
an alternate approach involving a variational (weak) formulation of the 
system equations along with estimates in the spirit of those found in 
numerous papers on finite element techniques in structural problems. In 
[10] such a treatment was given for damped cantilevered Euler-Bernoulli 
beams. In this presentation we outline this approach in the context of 
models for beams with tip masses and base acceleration. Full details of 
our results in this direction will be given in a more lengthy manuscript 
currently in preparation. 

Fundamental to our discussions is a conceptual framework in which 
one has a dynamical model with "states" u(t,x), 0 < t < T, x 6 ft, and 
"parameters" q(t,x), q 6 Q, where Q is an admissible class of 
parameter functions. The state system is an initial-boundary value 
problem involving a hybrid model (parameter dependent and coupled 
partial differential equations/ordinary differential equations). One is 
given observations (data) u ^ for u(t^,Xj) and seeks to solve the 
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optimization problem of finding parameters q" in the feasible parameter 
set Q which give a best (in the least squares sense) fit of the model 
to the data. 

We formulate this problem in an abstract setting with Hilbert state 

space V and parameter space Q. For computational purposes we then 

N M 

approximate V and Q by finite dimensional spaces V and Q 
respectively. We illustrate these ideas with a specific model and 

particular classes of approximations in our subsequent discussions here. 


II. The Identification Problem 

We consider a flexible beam of length £, spatially varying 
stiffness El and linear mass density p which is clamped at one end 
and free at the other with an attached tip mass of magnitude m (see 
Figure 2.1). 
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Using the Euler-Bernoulli theory to describe the transverse vibrations 
of the beam we obtain the partial differential equation (see [11], [12]) 

p(x)D^u(t,x) + D^EI(x)D^u(t,x) - D c(t,x)D u(t,x) + f(t,x), 

XX XX 

X € (0,£), t e (0,T), (2,1) 

for the transverse displacement u where a denotes the internal 
tension, f is the net externally applied transverse or lateral load, 
and = 3/3t,D x = 3/3x* Use of principles of elementary Newtonian 

mechanics (i.e., force and moment balance equations) yields the boundary 
conditions at the free end* From translational equilibrium we obtain 

mD^u(t,A) - D EI(£)D^u(t,£) = -a(t,A)D u(t,Jt) + g(t), t € (0,T), (2.2) 

where g is the net external force on the tip mass. In a similar 
manner, requiring rotational equilibrium, we have 

D^u(t,Jl) = 0, t e (0,T). (2.3) 

The geometric boundary conditions (zero displacement and zero slope) at 
the clamped end are given by 

u(t,0) =0, t € (0,T) , (2.4) 


and 
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D x u(t,0) =0, t € (0,T), (2.5) 

respectively. The initial conditions are in the form of initial 

displacement 

u(0,x) = <Kx), x € [0,£], (2.6) 

and initial velocity 

D t u(0,x) = t|>(x), x e [0,£], (2.7) 

In order to characterize solutions to the hybrid system (2.1) - 

(2.7) of ordinary and partial differential equations, boundary and 
initial conditions, we formally represent it as an abstract second-order 
system. Consider 


^0 D t U (t) + Ag U (t) 

= B Q (t)G(t) + F(t) , 

t € (0,T) , 

(2.8) 


A 

A 


YgU(t) = 0 

Yju(t) = 0 

Y 2 u(t) = 0 

(2.9) 

at x = 0, 

at x = 0, 

A A A 

at x = £, 

A 



u(0) = ♦, D t u(0) 

= 

(2.10) 


where u(t) = (u(t,£), u(t,»)), and the operators Mq, Aq, 8 Q (t) 
Y^> i = 0,1,2 are defined by (for v = (v(Jt),v) € R x H®(0,i)) 


and 
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MqV » (.mv(£), pvj, 

A q v = (.-D x EI(Jl)D^v(iO, D 2 EID 2 v} 

B 0 (t)v = (-a(t,£)D x v(A), D x oD x vJ 

and 

Y i V = D x V > 1 = °» 1 » 2 » 

F(t) = (g(t),f(t,*)J, <l> = (.<(>(£), 4>J, and • 

Define the Hilbert space H by 

H = R x H°(0,JO 

with innerproduct 

<(n,<|>),(5>'lO> H = ns + <<M> 0 

where {H* - , <• ,*>^) denote the usual Sobolev spaces together with the 
usual Sobolev inner products. Let V be the Hilbert space defined by 

V = {(n,4>) € H : * € H 2 (0,i), <j>(0) - D+(0) - 0, n - ♦(*)} 

together with the innerproduct 
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<i,i> v = <D 2 «j> , d 2 *> 0 , 

A A 

where <t> = (.<KJt),<t>J > t = (.’l>(£)»'l'J € V. It is easily shown that V is 
dense in H and choosing H as our pivot space we have the continuous 
embeddings V c H c V' where V' is the space of continuous linear 
functionals on V. 

Of particular interest to us here will be the notion of a weak 
solution to (2.8) - (2.10). Interpreting the derivatives in the 
definitions of the operators Aq, 8q( 0 and Y^ in the distributional 
sense, we rewrite (2.8) — (2.10) in variational form as 

<M 0 D 2 u(t),e> H + a(u(t) ,0) = b(t)(G(t),0) + <F(t) , 0> H * (2.11) 

0 € V, t € (0,T) 

u(0) = {, D t G(0) = } (2.12) 

where the sesquilinear forms a and b(t) on V x V are defined by 

a(i,i) = <EID 2 <(»,D 2 t> 0 

and 

b( t )(i,i) = -<oD x «t.,D x tj)> 0 

respectively, and the H innerproduct is interpreted as the duality 
pairing between V" and V (see [13], [14])) wherever appropriate. 
Under the assumption that El, p CL^O.Jl), a C [0,T] ,H*(0,£) J , 


f 6 L 2 ([ 0 ,T] > H®(0,£.)) and g € 1 < 2 ( 0 ,T), it is not difficult to show 
(see [15]) that the system (2.11), (2.12) admits a unique solution u 
with values in V and which satisfies u € C([0,T],v), D fc u € c([0,T],h) 
and Dj.u 6 [0,T] ,V'J . The existence of strong solutions can be 

demonstrated by rewriting ( 2 . 8 ) - ( 2 . 10 ) as an equivalent first-order 
system and using linear semigroup theory [16], [17] and evolution 
operators. The details involve standard ideas for evolution systems 
such as those found in [17] and under additional regularity assumptions 
(e.g., El € H 2 (0 ,£ ) , a €C 1 ([0,T], 1^(0,*)), 

$ € V fl {(v(£),v) | v € H 4 (0,£), D 2 v(£) = 0} , t e v), one can argue 
existence of strong solutions (i.e., u 6 c([0,T],v), D t u € c([0,T],h), 
Djru 6 ^([OjTJjH) with u satisfying (2.8), (2.9) almost everywhere 

on [0,T]) with sufficient smoothness to carry out the convergence 
arguments underlying the results presented in the next section. 

In formulating the identification problem, for ease of exposition 
we assume that we wish to identify the parameters m, El, p and a 
only. We do note, however, that our general approach is in fact 
applicable to a wider class of problems involving the estimation of the 
forcing terras and initial conditions (see [18], [19]). Let Q be a 

compact subset of £= R x 1^,(0, l) x 1^(0,*) x L 2 ( [0,T] ,H°(0,£)) . 

We assume that we have been provided with displacement observations 
{uXt^Xj) : 1=1, j=l , • • • ,v} at times € [0,T] and positions 
Xj € [ 0 ,£] and formulate the identification problem as a least squares 


fit to data: 
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(ID) : Find q = (m,EI,p,o) € Q which minimizes 

a U v _ 2 

J(q;u) = l l |u(t.,x.) - u(t.,x )| , 
i=l j=l J J 

A 

subject to u(t) = (u(t,£)>u(t,*)) being the solution to (2.11), (2.12) 
corresponding to q € Q. 


III. The Approximation Scheme 

Our approximation scheme is based upon the use of a standard 
(finite element) Galerkin approach to construct a sequence of finite 
dimensional approximating identification problems. For each 
N = 1,2, # »*, let c V be a finite dimensional subspace of H. 

Let denote the orthogonal projection of H onto with respect 

to the H innerproduct. The Galerkin equations for the system (2.11), 
(2.12) are 

<M 0 D^u N (t),e N > H + a(G N (t),e N ) = b(t)(G N (t),e N ) + <F(t),e N > H , (3.i) 

e N e v N , t e (o,t), 

G N (0) = P N J, D fc G N (0) = P N J>, (3.2) 

where u N (t) = (u N (t,i), u N (t,»)) € V N . The approximating 

identification problems then take the form 
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(IDN) : Find = (m^,EI^,p^,a^) € Q which minimizes J(q;u^) 

subject to u being the solution to (3.1), (3.2) corresponding to 

q € Q* 

Of particular interest to us here is a scheme involving the use of cubic 
spline functions. Let denote the (modified) cubic B-splines 

on the interval [0,£] corresponding to the uniform partition 
{0, ^ which satisfy Bj(0) = DBj(O) = 0, j = 1,2,***,N+1. 

Let V N = span{(Bja), B^)}^}. Then V N c v and (3.1), (3.2) take the 
form 

+ A N u N (t) - B N (t)w N (t) + F^t), e € (0,1) 

w N (0) - [W^r 1 ; N (o> - iwV 1 

where 

:■<») - 1 1 At>(»?(i>.*;). 

j _1 J J J 

= dbJoOBjW) + / pB^ bJ, 

[A N ] ±j = / EID 2 B^ D 2 Bj, 

[B N (t)J = / o(t,*)DB^ DB N , 

1J 0 J 

£ 

[F N (t) 1 = g(t)B^Ol) + / f(t,») B ?» 

x o 1 


(3.3) 

(3.4) 
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[Wp] i = <KjObJ(jO + / 


tw^] = *(i)B?U) + / *B?, 

1 o 


Iw \j 


= bJ(DbJw) + J b^ Bj n , i,j 


Our convergence results for the cubic spline approximation schemes are 
summarized in the following two theorems. 


Theorem 1. Suppose {q^} c Q with q _as^ N + <*>. Suppose 

further that u(q), the solution to (2.11), (2.12) corresponding to 



Theorem 1 can be established using approximation properties of 
cubic splines (see [20]) and variational arguments which are similar in 
spirit to those found in [21] for second-order hyperbolic systems and in 
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[10] for damped cantilevered beams without tip mass. Continuous 
dependence of J and u and compactness of Q allow us to conclude 
that problem (IDN) admits a solution. The existence of a convergent 
subsequence also follows from the compactness of Q. Finally an 
application of Theorem 1 yields 

_ - _ _N V „N _N, »N, 

j( q»u(q)) = lim j(q ,u (q )) < lim j(q,u Aq)) = j(q,u(q)) 
k->“ k-*- 00 

for all q € Q and Theorem 2 is thus proven. 

Although the state equation in problem (IDN) is finite dimensional, 
the admissible parameter space Q is a function space and hence the 
minimization of J is over an infinite dimensional space* We briefly 
indicate a means of overcoming this difficulty which involves the 
introduction of a second level of approximation into our scheme. A 
detailed discussion of these ideas along with several numerical examples 
for problems with parabolic, hyperbolic, and simple Euler-Bernoulli 
equations can be found in [22], [23], and [7] respectively. 

For each M = 1,2,»»*, define the set Q M c £ by Q M = I M (Q) 
where I** is a mapping which satisfies 

M 

Ql. I : Q + £ is continuous 

Q2. I M (q) + q as M ^ 00 uniformly in q for all q 6 Q. 

The approximating identification problems now take the form 
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(IDNM) : Find = (m^,EI^,p^,oj^) € which minimizes J(q;u N ) 

subject to u N being the solution to (3.1), (3.2) corresponding to q. 

Typically, the spaces and the mappings 1^ are realized using 

finite dimensional spaces of interpolating linear or cubic spline 

functions. In this case under sufficient regularity assumptions on Q, 

it can be shown that conditions Q1 and Q2 above are satisfied. 

Using conditions Ql, Q2, and the compactness of Q, one can readily 

establish a convergence result analogous to that given in Theorem 2. 

Specifically, if {qjJ} is any sequence of solutions to the problems 

r _ N k. _N k _ 

(IDNM), there exists a convergent subsequence i^m.I with q^ q 

as N v “, M. +• °°, where q is a solution to (ID). 

K J 


IV. A Numerical Example 

We present a representative example to illustrate some of the 
numerical results we have obtained using the methods outlined above. 

Further details and other numerical findings will be presented 

elsewhere. 

We consider the problem of estimating the spatially invariant 

stiffness El and linear mass density p of a cantilevered beam of 

length & = 1 with an attached tip mass at the free end of unknown 
magnitude m which is also to be identified. We also assume that the 
entire system is subjected to a time varying base acceleration ag(t). 
The internal tension o is then given by (see [9], [12]) 


T 
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o(t,x) = - a Q (t)(p(£-x) + m). 

The system was assumed to be initially at rest C <t> “ = 0) and then 

acted upon by the distributed transverse load 

f(t,x) = e X sin2irt 

and point load at the tip 

g(t) = 2e fc . 

The base acceleration aQ was taken to be 

! 1 0 < t < 1.5 

0 t > 1.5 

"Observations" (i.e., displacement values to be used as data in the 
inverse algorithm) at positions Xj = .75, .875, 1.0 at times tj = .5, 
1.0, ••*,5.0 were generated using the "true" values of the parameters 
m = 1.5, El = 1.0, and p = 3.0, the first two natural modes of the 
unforced, unaccelerated system and a standard Galerkin scheme. The 
approximating optimization problems were solved using a Levenberg- 
Marquardt iterative steepest descent method. "Start up" values for the 
parameters to be estimated were chosen as mg = 1.7, EIq * .7 and 
Pq = 2.7. The initial value problem (3.3), (3.4) was solved at each 
iteration using a variable step size Adams predictor corrector method. 
The system did not appear to be stiff. Our results are summarized in 


Table 4.1 below. 
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Table 4.1 


N 

El* 

-N 

P 

-N 

m 

J* 

CPU (min/sec) 

2 

1.0016 

3.0997 

1.4793 

.17 x 10" 4 

0/18.01 

3 

1.0012 

3.0636 

1.4873 

1 

o 
1— 1 

X 

00 

i-H 

• 

0/35.93 

4 

1.0009 

3.0414 

1.4921 

.19 X 10" 4 

1/26.19 

5 

1.0006 

3.0306 

1.4944 

.19 x 10“ 4 

4/15.87 

6 

1.0012 

3.0319 

1.4944 

.35 x 10" 4 

5/21.34 
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